Critical Metal Phase at the Anderson Metal-Insulator Transition with Kondo Impurities 
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It is well-known that magnetic impurities can change the symmetry class of disordered metallic systems by 
breaking spin and time-reversal symmetry. At low temperature these symmetries can be restored by Kondo 
screening. It is also known that at the Anderson metal-insulator transition, wave functions develop multifractal 
fluctuations with power law correlations. Here, we consider the interplay of these two effects. We show that 
multifractal correlations open local pseudogaps at the Fermi energy at some random positions in space. When 
dilute magnetic impurities are at these locations, Kondo screening is strongly suppressed. We find that when 
the exchange coupling J is smaller than a certain value J*, the metal-insulator transition point extends to a 
critical region in the disorder strength parameter and to a band of critical states. The width of this critical region 
increases with a power of the concentration of magnetic impurities. 



Fifty years after its proposal, the Anderson metal-insulator 
transition (AMIT) of disordered noninteracting electrons ^ 
continues to be intensively studied 12. S H ISfl ■ The AMIT is 
a quantum phase transition of second order where the local- 
ization length diverges with a critical exponent v. The critical 
state is multifractal and characterized by a wide distribution of 
wave function amplitudes with a log-normal shape jH. One 
aspect which remains less understood is the interplay between 
the AMIT of conduction electrons and the dynamics of the 
spin of magnetic impurities. These can be induced by energy 
levels below the Fermi energy such as the d-levels of transition 
metal impurities |6] or by localized electronic states such as 
the dopant levels in semiconductors f^. Magnetic moments 
can enhance the spin susceptibility and the specific heat as 
has been observed in Si:P close to the AMIT When lo- 

cal magnetic moments interact with the conduction electrons 
by an exchange interaction, they can break both spin degener- 
acy and time-reversal symmetry (TRS) of the conduction elec- 
trons. Therefore, they are expected to change the symmetry 
class of the electronic system from orthogonal to unitary jioll . 
changing the critical exponent v, the critical electron density 
ric, and the critical disorder strength Wc at which the tran- 
sition occurs im \vh . At sufficiently low temperatures, an 
additional effect comes into play. The antifeiTomagnetic ex- 
change interaction between spin- 1/2 local moments and the 
conduction electrons in a metal leads to a correlated electron 
state where a Kondo singlet is formed, screening the local mo- 
ments at temperatures below the Kondo temperature Tk- In 
this limit, the magnetic susceptibility x approaches a constant 
value, as indicated in the left inset of Fig. [T] and the magnetic 
moments cease to break the TRS of the conduction electrons. 
The situation is quite different in the insulating phase, where 
local spectral gaps A/ prevent the development of the Kondo 
screening whenever the exchange coupling J between local 
moments and conduction electrons is below a critical value 
1 13]. In this case, the unscreened, free magnetic moments 
(FMMs) yield a low-temperature susceptibility x ^ 1/T, as 
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FIG. 1 : (Color online) The fraction of free magnetic moments Pfm 
[Eq. ( I13b 1 in a three-dimensional disordered electron system as func- 
tion of the exchange coupling J (in units of band width D) and the 
disorder strength W (in units of the critical value W^). Dashed line: 
J^(W^), Eq. the critical exchange coupling caused by local spec- 
tral gaps due to Anderson localization. For J < J* ~ 0.3Z), there 
is a critical metal region for W° < < M^c( J), Eq. d. Insets: 
temperature dependence of the local magnetic susceptibility x of ^ 
single magnetic impurity in the metallic and in the insulating phase, 
showing Kondo screening and free moment behavior, respectively. 



indicated in the right inset of Fig. [T] As a result, TRS may be 
partially broken in the insulating phase. The symmetry class 
at the AMIT is determined by the subtle interplay between the 
Kondo screening and the Anderson localization. 

In this Letter, we aim at determining the density of FMMs 
in the vicinity of the AMIT in the low-temperature limit 
T <C {Tk) for a system of dilute magnetic impurities with an 
average Kondo temperature {Tk). An essential aspect of the 
problem are the multifractal fluctuations of the eigenfunctions 
at the AMIT. They lead to the local suppression of wave func- 
tion amplitudes at random locations 114]- Multifractal states 
are power-law coiTelated over a large energy interval Ec of 
the order of the band width D [ISL Il6|] . We show that these 



spectral coiTelations induce local pseudogaps which, in turn, 
prevent the Kondo screening of a sizeable fraction of local 
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moments even when the exchange couphng exceeds (see 
Fig. [U. When the exchange couphng J is sufficiently weak, 
these FMMs can break TRS and shift the transition to stronger 
disorder amphtudes, Wc{J), as well as lower the critical ex- 
ponents on the insulating side to their unitary values. Com- 
ing from the metallic side, the transition occurs at the smaller 
disorder strength since the magnetic impurities become 
Kondo screened. We conclude that the critical point is ex- 
tended to a critical region, as shown in Fig. [T] 

The most general model for the formation of magnetic im- 
purities is the Anderson model [6] describing a localized level 
with energy and on-site Coulomb repulsion U hybridizing 
with electrons in the conduction band of a Hamiltonian Hq 
which may include the random potential caused by nonmag- 
netic impurities ^(r). Using as basis the eigenstates V'n with 
energy £"„ of Hq and the density operator fin^^, the Anderson 
Hamiltonian can be written as 



Ha ^ 



nn.a + Erf 



fld.a + U 



nd,+nd, 



J2{tndC+,d^ + n.c.) 



(1) 



Here, fid^a is the density operator in the impurity level. 
The hybridization amplitude tnd is proportional to the eigen- 
function amplitude at the position of the magnetic impurity 
tpni^) ^nd the localized orbital amplitude (pdiO): tnd = 
til'*(0) (j)diO). One can employ a Schrieffer- Wolff transfor- 



mation lllTl 11811 formulated in terms of eigenstates tpn 1.19.1 to 
generate an s-d Hamiltonian with exchange couplings 



J n 



^nd ^n'd 
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U + Ed — En' En — (d 



(2) 



and an additional potential scattering term with amplitude 
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U + ed — En' En — ed 



(3) 



Note that Knn' vanishes for all n, n' when \En.n' ~ Ep\ <C U 
and = Ep — C//2, where Ep is the Fermi energy. For ar- 
bitrary wave functions with small amplitude at the position 
of the magnetic impurity are hardly modified by this potential 
scattering term since Knn' ^ '0n(O)V'n'(O) we can disre- 
gard it when deriving the density of FMMs. Thus, we retain 
only the exchange couplings Jnn' = Jipn{^)i'ri'{Q)- 

Critical exchange coupling. - Recently, it was shown that 
the distribution of the exchange couplings J„„/ results in a 
wide bimodal distribution of the Kondo temperature and free 
magnetic moments 12111 . This was at first obtained 

by solving the 1-loop Nagaoka-Suhl equation (NSE) [I22I] of 
the Kondo model in the representation of the eigenstates of 
Hq. It was confirmed with the numerical renormalization 
grou p and the continuous time quantum Monte Carlo method 
112 111 . These nonperturbative methods show that the NSE un- 
derestimates the density of FMMs. Building on that work, we 
can get the lower limit of the critical exchange coupling Jc by 



setting T/f ^ in the NSE: 
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2N^^ \En'~Ep\ 

71—1 ' ' 



(4) 



where L denotes the linear size of the system, N ^ L'^ the 
number of states in the band, and d is the spatial dimension. 
For a clean metallic system with a flat band, the critical ex- 
change coupling vanishes logarithmically with the number of 
states N as Jc ^ D/ \nN, leaving no FMMs for any fi- 
nite value of J. On the other hand, if the eigenstates at the 
Fermi energy become localized with a localization length ^, 
finite local gaps of order A/ — {pS,'^)^^ appear at the Fermi 
energy, cutting off the Kondo renormalization 11 1311 (p is the 
average density of states). Thus, there are FMMs whenever 
A/ > Tk, or, equivalendy, J D/lnNj, where 

A^i — D/Ai is the number of localized states with a finite 
wave function amplitude at the magnetic impurity site |21]. 

Critical exchange coupling at the AMIT. - For d > 2 an 
AMIT exists and the localization length diverges at the mobil- 
ity edge £;me as a power law, ^{E) ^ {\E - EmeI/Ec)'". 
Another way to probe the transition region is to vary the dis- 
order strength W. At a critical strength Wc, the localization 
length diverges with the same exponent i^. Thus, one can draw 
a quantum phase diagram with J versus W where the EMM 
phase is limited by a line defined by the function 

J^(W)/D = [lydlniEJlW - Wei)]"' , (5) 

as shown by the dashed line in Fig. [T] 

Local Pseudo Gaps. - Equation ^ does not take into ac- 
count the multifractality and the critical correlations between 
wave functions at different energies at the AMIT |15]. The 
amplitude of multifractal states is suppressed at some random 
positions below their typical value, scaling as L^" with a > d 
(i.e., decaying faster than extended states). Correlations be- 
tween wave functions at different energies can then open wide 
local pseudogaps. These correlations can be quantified by 
spatially integrating the correlation function of eigenfunction 
probabilities associated with two energy levels distant in en- 
ergy by uinrn ^ En - E„i, namely [16], 



Cn 



= L" 



< A, 



{Ej\uJnrn\)\ A < \uJn,n\ < E^, (6) 
{Ec/\uJnrn\)'^, \l^nm\ > Ec, 

where < 7 < 1 . For LOnm < Ec the value of Cnm is en- 
hanced in comparison to the plane-wave limit Cnm — 1- The 
distribution function of a multifractal wave function is radi- 
cally different from the Porter- Thomas distribution of metallic 
states since the moments of eigenfunction intensities |7/';(r)p 
obey the power-law scaling Pq = L'^{\tpi{r)\'^'^) i"'^', 
where Tq is the multifractal exponent of the q-th moment. 
The corresponding distribution function is known to be log- 
normal, in good approximation [5], 

P{a)^L 53^, (7) 
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where a = — \n\tpi{r)\'^ / \n L. ao > d is related to the mul- 
tifractal exponents, Tq = d{q — 1) + (ao ~ d)q{l — q), and to 
7 — 2{ao — d)/d. In order to find we need the joint distri- 
bution of two eigenfunction intensities |-0i(r)P and |?/'„(r)p, 
I 7^ n which is also log-normal, 

with f{a) = d — {a — oq)'^ /2d'y. Imposing that Eq. ([8]l 
is consistent with Eqs. (|6]l, we determine the coefficients as 

ttki = (1 + Vl + 46fc;)/2 and bki = gki/igl, - 1), where 



1, 



Jkl 



_ lii{\ujik\/Ec) , . ^„ 

" d\nL 1 2/7, Ic^/fcl >i?e. ^ ^ 

The conditional probability to have a; when q;„ = a is then 

[a; -a0+3i„(o 



P{a) 



(10) 



Averaging Eq. (|4| with the conditional probability Eq. ( fTOl ), 
we find that, at a site where = a with ~ -Bf. the 
critical exchange coupling Jc is determined by the equation 



El —Ef -^Ec 



{\Ei-Ef\/E, 



NlEi - Ef 



(11) 



The factor {\Ei — Ef\/Ec) ^ corresponds to a local pseu- 
dogap with exponent (a — ao) /d at positions where a > ao- 
At such sites the critical exchange coupling is proportional to 
that exponent, in the continuum limit A — > / dEi, 



Jcja) 
D 



a — ao 



1 — exp 



(12) 



Notice that when a takes its typical value ao, we get Jc(ao) — 
y^2dj/d\n L, which decreases slowly with increasing L. At 
sites where a < ao, the LDOS is larger than its typical 
value and Jc vanishes. The exchange coupling ./i'^^ at which 
there is on average only one EMM in the whole system is 



X 



(1) 



2jD, where a+ is obtained from the 



condition P(a+) = 1/L , yielding a+ 
r(2 



ao 



2-fd. Sim- 
ilarly, at J^^' = ^D, one finds 0{\/ N) sites with FMMs. 
A finite density of FMMs at the AMIT is only found at very 
small exchange couplings, ,/ < Jc(ao), and these values are 
vanishingly small for L ^ oo. 

Pseudogaps and FMMs in the insulating phase. -While we 
conclude that there is no finite density of FMMs at the AMIT 
in the infinite-volume limit, the situation changes on the insu- 
lating side of the transition. First, a hard spectral gap of order 
A/ quenches all magnetic moments below an exchange cou- 
pling Jjf , given by Eq. (|5]l. However, there are still power-law 
spectral correlations between wave functions located within 
the same localization volume One can take these correla- 
tions into account by noting that the wave function intensities 



within a localization volume still have a log-normal distribu- 
tion with a — > a^ = — In It/i^/ In^. For the evaluation of 
Jc we therefore should substitute L by the localization length 
{{W) ^ (W - WcY". Thus, for fixed J, the density of 
FMMs depends on the localisation length as 



(13) 



Pfm = n-pyilnu = Erfc W-j^t^ 




Close to the transition, is large and Eq. ( fT3] l simplifies to 
Pfm ~ (W^ " W^c)''^'^\ decaying with an exponent k( J) — 
(vdl2^){JlD)\ 

Metallic phase. - On the metallic side of the AMIT, the 
correlation length limits the range of multifractal corre- 
lations. One can imagine the extended states close to the 
transition as patches of multifractal states connected to each 
other by tunneling |23]. This yields the scaling IV'P'^ ~ 
{icIEf'^ic'^^'^'' - for L > ic. leaving no finite density 
of FMMs in the metalHc phase (W < Wc). 

Symmetry Class of the AMIT. - The magnetic scattering by 
FMMs changes the symmetry class of the conduction elec- 
trons from orthogonal to unitary when /Vts > 1, where 
T) is the diffusion constant, 1/ts is the magnetic scatter- 
ing rate, and ^ denotes the correlation length on the metal- 
lic side, or the localization length on the insulating side of 
the transition, respectively ifiol Til , 24]. The magnetic scat- 



tering rate at zero temperature is bounded from below by 
27rni?M'S'(5'+l)(J/£')^p(ei?,r)/p2. Coming fromthe metal- 
lic side of the transition, npM vanishes due to the Kondo 
screening and the AMIT occurs at the orthogonal critical value 
for time-reversal symmetric systems, W^. On the insulating 
side of the AMIT, the density of FMMs is finite, as given by 
Eq. ((ni). The LDOS at the position of the FMMs scales as 
p(eF,r) - ^d-a(.j)^ where, a(J) = ao + dJJD. There- 
fore, 1/ts is finite and, following Refs. ifTH 24], we can 
get the critical disorder amplitude Wc{J) from the condition 



DT,ij) = ^{Wc{j)~Wc''y 



WciJ)~Wc^ 



2'^{TrWciJ)J)^nMS{S + l) 
dD^ 



niJ) 



(14) 

where l/ri{J) ~ 2v — {dv/2^){'-^ + J / D)"^ . For J exceeding 
J* = {2yJ^/d — "f)D, the critical disorder amplitude Wc{J) 
approaches its orthogonal value W^. However, for smaller 
exchange couplings, J < ,/*, a paradoxical situation appears: 
The Critical Metal Phase. - The position of the critical 
point Wc depends on the direction from which the AMIT is 
approached: Coming from the metallic side, the Kondo effect 
screens all magnetic moments at low temperature T ^ Tk, 
so that TRS is not broken and the AMIT occurs at W^. Com- 
ing from the insulating side, there is a finite density of FMMs 
which break TRS and lead to an AMIT at the stronger disor- 
der amplitude Wc( J) given by Eq. (fT4l i. Consistency requires 
then that these two distinguished critical points are connected 
for intermediate disorder strengths < W < Wc{J) by 
a critical metal region. Accordingly, the mobility edge is 
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extended to a critical band whose width is a function of J 
and um- The resulting zero-temperature quantum phase dia- 
gram is shown in Fig. [Tlwith the critical metal phase appear- 
ing below the tricritical point J*. The analytical expression 
for the density of FMMs compares quantitatively well with 
results obtained from numerical simulations of critical ran- 
dom banded matrices ll25[l . and are in qualitative agreement 
with numerical simulations of the three-dimensional AMIT 
lEoll . Furthermore, we have derived [25] the distribution of the 
Kondo temperature analytically taking into account multifrac- 
tality and critical correlations of the critical wave functions 
and found a wide bimodal distribution at the AMIT. In two di- 
mensions, the distribution is also bimodal with a finite width 
which increases with the second power of the disorder am- 
plitude, in qualitative agreement with numerical simulations 

lis El. 

A wide distribution of Kondo temperatures was argued pre- 
viously to cause a finite electron dephasing rate at low tem- 
peratures 1 19] and to lead to deviations from the Fermi-liquid 
behavior [7,^- Here we consider the zero-temperature limit 
and find that on the metallic side of the AMIT the density of 
FMMs vanishes due to the Kondo screening. However, a finite 
density of FMMs survives on the insulating side of the AMIT. 
A new critical metal phase arises due to the interplay between 
the Kondo effect and the multifractality of the wave functions. 
One experimental consequence is the divergence of the dielec- 
tric susceptibility, x(T = 0) - ^ (WdJ) - W)-^" at 
the disorder amplitude Wc{J) (or at a critical electron den- 
sity nc{J)), while the zero temperature limit of the resistivity 
diverges with the correlation length as p{T = 0) ^ ^ 
{W — W^)~'^ at the weaker disorder amplitude (or at 
a larger critical density , accordingly). This difference 
Wc{J) — increases with the concentration of Kondo im- 
purities as Eq. ( fT4] i. maximally by about 10%fl2']. We ex- 
pect the critical metallic phase to be observable in materials 
where both the AMIT and the Kondo effect are present simul- 
taneously at experimentally accessible temperatures, such as 
in amorphous metal-semiconductor alloys [3, 26] with dilute 
magnetic impurities ll27ll . or in doped semiconductors, such as 
Si:P where thermopower measurments are consistent with the 
presence of Kondo impurities with (Tk) ~ li^LlSJ- 
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